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Abstract.  A refinement  to  the  theory  published  by  Finn  et  al.  is  presented  here.  Compression  effects 
are  taken  into  account  by  a rigorous  definition  of  the  plasma  length  and  by  modifying  the  expression 
of  the  velocity  field.  The  perturbation  of  the  plasma  length  is  calculated  exactly  by  a suitable  Green 
function.  Growth  rates  and  real  frequencies  of  the  unstable  mo  = 1 mode  are  compared  with  the 
experimental  values,  showing  a good  agreement  when  compression  effects  are  strong  (i.e.,  for  short 
traps). 


INTRODUCTION 

The  niQ  = 1 diocotron  instability  is  one  of  the  open  problems  in  non-neutral  plasma 
theory.  According  to  the  classical  theory,  the  m e = 1 diocotron  spectrum  does  not  present 
discrete  unstable  modes  (regardless  of  the  equilibrium  density  profile)  [1],  while  the 
continuum  spectrum  can  only  produce  an  algebraic  growth  proportional  to  tx! 2 [2].  On 
the  contrary,  experiments  show  that  the  linear  growth  of  the  mode  is  exponential  [3, 4]. 

Finn  et  al.  [5]  provided  a justification  of  the  m$  = 1 diocotron  instability  in  terms 
of  compression  of  the  plasma  column  due  to  the  confining  potential  along  the  axial 
direction  of  the  trap.  However,  the  theory  by  Finn  et  al.  neglects  the  finite  electron 
temperature  and  the  effects  of  the  perturbation  of  the  plasma  length  in  the  axial  direction 
are  treated  in  a simplified  manner.  These  two  issues  are  treated  rigorously  in  a recent 
work  by  the  Authors  [6]. 

In  particular,  the  new  model  differs  from  the  one  in  Ref.  [5]  for  three  main  aspects: 

1.  a rigorous  definition  of  the  length  of  the  plasma  column  is  used; 

2.  an  effective  electrostatic  potential  is  introduced,  in  order  to  calculate  the  E x B drift 
on  a string  of  variable  axial  density  (finite  temperature  effects  are  included); 

3.  the  perturbation  of  the  plasma  length,  induced  by  density  variations,  is  calculated 
in  terms  of  a suitable  Green  function. 

In  the  present  work  the  theory  presented  in  Ref.  [6]  is  reviewed  and  comparisons  with 
experimental  data  provided  by  Kabantsev  and  Driscoll  [4]  are  presented. 
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THE  PHYSICAL  MODEL 


A Penning  trap  confining  an  election  plasma  is  considered  in  the  following.  The  central 
electrode,  which  extends  between  z = —Lc/2  and  z — Lc/2,  is  grounded,  while  the  end 
electrodes  (each  having  length  Ls ) are  at  a negative  potential,  —V.  The  radius  of  the  trap 
is  indicated  by  /?w  and  the  gap  between  the  central  and  the  end  electrodes  is  Lg. 

The  main  difference  of  the  present  model  with  respect  to  the  classical  one  is  regarding 
electrons  as  string  of  variable  length  and,  consequently,  of  variable  density  (while  in  the 
classical  approach  they  are  considered  as  strings  of  uniform  density).  Since  the  plasma 
length  changes  radially  (due  to  the  confining  potential),  the  columns  are  compressed  or 
rarefied  as  they  move  radially.  This  phenomenon  is  known  as  compression  effect. 

The  model  can  be  deduced  starling  from  the  following  two  hypotheses: 

1.  Due  to  the  presence  of  an  axial  magnetic  field,  in  the  characteristic  time  scale  of 
the  diocotron  instability  the  Maxwell-Boltzmann  distribution  is  reached  only  along 
the  z direction  [7]  and  the  particle  density  is 


n (r,Q,z,t) 
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where  © = kT /e  is  the  electron  temperature  expressed  in  Volts,  while  nc  ( r,0,t ) and 
(j)f  (/-, 0, /)  are  respectively  the  electron  density  and  the  electrostatic  potential  in  the 
central  region  of  the  trap,  i.e.,  for  |z|  < 0 /2. 

2.  Electrons  move  in  the  transverse  plane  due  to  the  E x B drift  velocity: 


vj.(r,0,z,r) 
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The  continuity  equation 

Instead  of  nc  (r,  0,t),  experiments  yield  the  line  integrated  density,  o(r,  0,/),  given  by 

a(r,Q,t)=  f n(r,Q,z,t)dz  = nc(r,Q,t)L(r,Q,t),  (3) 

J -oo 

where  the  effective  length  of  the  plasma  column,  L(r,Q,t),  has  been  defined  as 
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Averaging  the  continuity  equation  for  n ( r,Q,z,t ) over  z,  one  obtains  a new  continuity 
equation  for  the  line  integrated  density 


^ + Vx-(VJo)  = 0 

with  the  velocity  given  by  a suitable  z-averaged  E x B drift. 


(5) 
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The  z-averaged  E x B drift 


The  present  theory  differs  from  the  one  by  Finn  et  al.  [5]  in  the  evaluation  of  Vx- 
In  fact,  the  velocity  is  here  calculated  in  terms  of  the  E x B drift  estimated  on  a whole 
string  (a  similar  assumption  has  been  employed  by  Fine  and  Driscoll  [8])  and  it  is  given 

by 
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Equation  (6)  can  be  rewritten  as 


V±  = ^-ez  x Vx4>eff, 

a 0 

where  the  effective  potential,  <|)eff(r,0,f),  is  defined  as 
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<t>eff  = <l>6-  + ©log 
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The  model 


The  model  is  completed  by  the  Poisson  equation  for  <J).  The  set  of  equations  can 
be  written  in  a more  compact  way  by  defining  the  quantity  'F(r,01z)f)  = <|>(r,0,z,f)  - 
(j>c  (r,  0, r),  which  represents  the  correction  to  the  potential  in  a trap  of  finite  length,  with 
respect  to  the  one  in  a trap  of  infinite  length.  Finally,  the  following  system  of  equations 
is  to  be  solved: 


dt 
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TABLE  1.  Confining  potential  and  geomet- 
rical parameters  of  the  experiments  performed 
by  Kabantsev  and  Driscoll  1 4] . 


V\V] 

L,  [cm] 

Lf,  [cm] 

Ls  [cm] 

A 

36 

34.6 

0.175 

6.8 

B 

80 

20.7 

0.175 

13.8 

C 

too 

7 

0.175 

6.8 

FIGURE  1.  Experimental  line  integrated  density  from  experiments  by  Kabantsev  and  Driscoll  1 4] . 


A modified  Rayleigh  criterion  for  linear  stability  and  the  conscrvalion  of  canonical 
angular  momentum  and  of  the  energy  of  the  system  can  be  deduced  from  system  (9) 
[6].  The  criterion  proves  that  when  ao(r)  is  a monotonic  function,  the  system  is  linearly 
stable. 


COMPARISON  WITH  EXPERIMENTAL  RESULTS 

Growth  rates  calculated  with  the  present  theory  have  been  compared  with  those  resulting 
from  recent  experiments  performed  by  Kabantsev  and  Driscoll  [4],  For  each  experiment 
the  electron  temperature  is  approximately  1 V and  the  axial  magnetic  field  is  Bo  — 0.4  T. 
Three  different  cases  (labelled  A - B-C)  have  been  considered.  Table  1 summarizes  the 
confining  voltage  and  the  geometrical  parameters  of  the  trap  for  these  cases. 

For  each  case,  the  initial  z-integrated  density,  Go(r),  is  reported  in  Fig.  1.  From  such 
measured  density  distributions,  the  corresponding  profiles  for  nc  ( r ) and  L(r)  have  been 
calculated,  by  solving  Eqs.  (9)  iteratively  for  the  initial  equilibrium.  The  results  of  the 
calculations  are  presented  in  Figs.  2 and  3.  Growth  rates  of  the  instability,  y,  and  real 
frequencies,  /,  for  the  mQ  = 1 mode  have  been  evaluated  by  linearizing  system  (9). 
Table  2 summarizes  the  results. 
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FIGURE  2.  Equilibrium  density  profile  calculated  from  experimental  Oq. 


A few  considerations  are  in  order.  It  can  be  noticed  that  in  a short  trap  (case  C)  com- 
pression effects  are  very  strong  and  the  equilibrium  length  profile  varies  considerably 
with  r.  On  the  contrary,  for  longer  traps  (cases  A and  B)  the  equilibrium  length  profile 
has  a small  variation  with  respect  to  its  value  for  r = 0 (10%  inside  the  plasma  radius). 
When  compression  effects  are  small  (i.e.,  for  long  traps),  the  growth  rates  calculated 
with  the  present  theory  are  not  comparable  with  the  corresponding  experimental  values. 
Viceversa,  when  compression  effects  are  strong  (i.e.,  for  short  traps)  the  agreement  be- 
tween theory  and  experiment  is  remarkable.  When  compression  effects  are  small,  the 
disagreement  between  theory  and  experiment  is  probably  to  be  explained  in  terms  of 
kinetic  effects,  as  pointed  out  in  Refs.  [9]  and  [10]. 
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TABLE  2.  Real  frequency  and  growth  rate  (in  kHz)  for 
the  unstable  »iq  = 1 mode 


/ 

experiment 

Y 

V// 

/ 

theory 

Y 

y/f 

A 

17.9 

6.0 

0.34 

17.0 

0.99 

0.058 

B 

23.2 

5.6 

0.24 

21.0 

0.93 

0.044 

C 

42.7 

10.6 

0.25 

36.1 

8.5 

0.24 
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